Orbital occupation, local spin and exchange interactions in V2O3 
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We present the results of an LDA and LDA+U band structure study of the monoclinic and the 
corundum phases of V2O3 and argue that the most prominent (spin |) models used to describe 
the semiconductor metal transition are not valid. Contrary to the generally accepted assumptions 
we find that the large on site Coulomb and exchange interactions result in a total local spin of 1 
rather than | and especially an orbital occupation which removes the orbital degeneracies and the 
freedom for orbital ordering. The calculated exchange interaction parameters lead to a magnetic 
structure consistent with experiment again without the need of orbital ordering. While the low- 
temperature monoclinic distortion of the corundum crystal structure produces a very small effect 
on electronic structure of V2O3, the change of magnetic order leads to drastic differences in band 
widths and band gaps. The low temperature monoclinic phase clearly favors the experimentally 
observed magnetic structure, but calculations for corundum crystal structure gave two consistent 
sets of exchange interaction parameters with nearly degenerate total energies suggesting a kind of 
frustration in the paramagnetic phase. These results strongly suggest that the phase transitions 
in V2O3 which is so often quoted as the example of a S=| Mott Hubbard system have a different 
origin. So back to the drawing board! 



The V2O3 system has been a topic of intense study 
for over more than fifty years by both theoreticians and 
experimentalists because of it's "rich" phase diagram. It 
undergoes a first order metal-insulator transition with a 
seven orders of the magnitude change in the electrical 
conductivity JO, which can be induced by temperature, 
pressure, alloying or by nonstoichiometry. It also exhibits 
an antifcrromagnetic inslulator to paramagnetic insula- 
tor transition which also is first order and a first order 
paramagnetic insulator to paramagnetic metal transition. 
Many theoretical models have been put forward to under- 
stand the electronic and magnetic behavior of this com- 
pound. Goodenough proposed a model involving both 
itinerant and localized 3d orbitals (||| to describe the 
the V2O3 electronic structure. Various models have been 
suggested and worked out in some detail using mainly the 
Mott-Hubbard picture || [|. In fact V2O3 is nowadays 
used as the best studied example of a Mott Hubbard sys- 
tem with a semiconductor to metal transition. To explain 
the peculiar antiferromagnetic order (Fig.|l|) in the low 
temperature insulating phase (AFI) with pairs of paral- 
lel spins coupled antiferromagnetically in the basal plane 
the intriguing idea of orbital ordering in the presumably 
doubly degenerate E g orbitals was suggested 0-^] and 
revived recently ]hH to explain new neutron scattering 
results Also within the context of infinite dimen- 

sion calculations to describe strongly correlated systems 
V2O3 is now used as a good example of the success of this 
approximation However, the character of the phase 
transition and the nature of the ground state is still very 
controversial. Recent photocmission and X-Ray absorp- 
tion results [|l3| strongly suggest that the AFI ground 



state should not be described as a spin | antiferromag- 
netic Mott-Hubbard insulator as assumed in the above 
theories and, related to this, that the d orbital occupa- 
tion is quite different from that conventionally assumed. 
In this letter we present the study of the electronic struc- 
ture of V2O3 in LDA @|]| and LDA+U approx- 
imations. In contrast to the standard LDA, in LDA+U 
method the influence of the on-site d-d Coulomb interac- 
tion is included in Hamiltonian, as an effective potential 
which is different for electron removal than for electron 
addition which has been shown to be crucial to describe 
strongly correlated materials. 




FIG. 1. the AF structure of the low temperature AFI phase 
of V2O3. The gray and filled circles correspond to spin-up and 
spin-down orientations of the local magnetic moments on V 
ions. The definitions of notations used for the superexchange 
interactions along the various paths are also shown. 

The crystal structure in the low temperature AFI 
phase is monoclinic |l8[ , above T c this changes to the 
corundum structure. (T c ~ 150K [[T9[). The calculated 
densities of states for both structures are very similar 
showing that the electronic structure itself is little influ- 
enced by the lattice distortions so that the strong change 
in properties must be a rather subtle effect with respect 



1 



to the electronic structure. In Figj^ the partial densities 
of states (DOS) obtained in the LDA calculation for the 
monoclinic crystal structure of V2O3 are shown. They 
are very similar to those found by Mattheis's [^0| for the 
corundum lattice. As the monoclinic distortion of the lat- 
tice is not very strong, we can plot DOS's assuming ap- 
proximate trigonal symmetry. The V ions are somewhat 
off center in a slightly trigonally distorted octahedron of 
O ions. This distortion causes the otherwise 3 fold degen- 
erate tig 3d orbitals to split into a non-degenerate A\ g 
and double degenerate E g levels. In this representation 
the Ai g orbital has 3z 2 — r 2 symmetry in a hexagonal co- 
ordinate system, i.e. with z axis along c direction (V-V 
pairs), and E g orbitals are directed more towards the V 
ions in the basal plane. 
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FIG. 2. Partial densities of state for 3d states of V obtained 
in the LDA calculation. Occupation numbers are given per 
one orbital (and both spins), i.e. each of the E g orbitals has 
occupation of 0.81 electrons. The Fermi level is at zero energy. 

At first glance the LDA picture looks very simple and 
one may conclude that the conventionally used ideas con- 
cerning the splitting of the A\ g orbitals into bonding 
and antibonding partners because of a hopping integral 
between the V ions in the pairs along the c axis with 
two electrons per pair in the bonding orbital are con- 
firmed. However upon closer examination a different pic- 
ture emerges which is further supported by the LDA+U 
calculations below as well as the recent experimental 
XAS results (T^|. First we see indeed a rather broad 
A\ g band with a total width of a little more than 2 eV 
and an about ^ eV narrower E g band both straddling 
the Fermi energy and resulting in a metallic state. The 
A\ g band, although it shows some structure that might 
be interpreted as a bonding-antibonding pseudo gap, ac- 
tually exhibits a strong peaking above Ef and only rel- 
atively little weight below Ef. This looks very different 
from a bonding-antibonding splitting which should have 
exhibited a more symmetric structure about the center 



of gravity if dimer hoping integrals were dominating the 
problem. We also note that the total A\ g band width is 
only 2.5 eV which is considerably smaller than the ex- 
pected value of the Hubbard U of about 3 eV including 
the screening due to the strongly bonding e g electrons 
pl[ or about 4-5 eV without this screening channel [ p2[ . 
Such a small band width would invalidate a molecular or- 
bital like approach also for the Ai g orbitals. So even these 
results already cast doubt on the validity of the most 
commonly used starting point with a strong bonding- 
antibonding splitting of the A\ g orbitals, with the E g 
orbitals in this gap ending up with two spin-antiparallel 
electrons per c-axis pair in the A\ g bonding state leaving 
only one electron in an assumed narrow doubly degener- 
ate E g band . It was this starting point which led to 
the now very much used one electron per site Hubbard 
model for V2O3. 

The total occupation number of the A\ g band is only 
0.5 electrons per V atom. Another very important point, 
which can be derived from the LDA result is that the E g - 
A\g actual splitting due to the trigonal crystal field is as 
large as ~0.4 eV with the E g band center of gravity lower 
than the Ai g band center of gravity. This large splitting 
makes the E g E g configuration for the localized AFI phase 
more favorable than the A\ g E g configuration removing 
the orbital degeneracy. We can clearly see at this point 
that if we switch on the on-site Coulomb interaction U in 
our LDA+U calculations the Ai g states will be pushed 
up above the Fermi level and the resulting ground state 
will be EgEg. This implies that the ground state (AFI 
phase) will not be degenerate leaving no place for any 
orbital ordering in the usual sense. 

The above LDA results give us rough ideas which 
should be further checked. The LDA+U calculations 
can provide us with some more details about the elec- 
tronic structure of V2O3. We performed LDA+U cal- 
culations for different magnetic structures, namely the 
"real" AF (Fig.0), a "simple" AF (all the nearest neigh- 
bors arc antifcrromagnctically coupled), a "layered" AF 
(all the nearest neighbors in basal plane are antiferromag- 
netically coupled and the neighbor along hexagonal c axis 
is coupled ferromagnetically) and FM structure. In Fig. || 
the LDA+U partial DOS are plotted. Here the results 
for U=2.8 (eV) and J=0.93 (eV) are presented. These 
values of U and J were calculated by Solovyev et. al, 
taking into account the screening of tig interactions by 
e g electrons plj . The striking point of the LDA+U re- 
sult is that the electronic structure strongly depends on 
the magnetic structure. The band gap for instance in 
the real AF magnetic structure is 0.6 (eV) close to the 
experimental value [Fig.|(a)], while in the FM structure 
we find a half-metal[Fig.|(c)]. The nature of the DOS 
changes strongly even far below and and above Fermi 
level indicating the very strong sensitivity of the elec- 
tronic structure to the spin structure. However, the oc- 
cupation numbers for all magnetic structures are ~ 0.95 



2 



for each of E g orbital, and ~ 0.25 for A\ g , which for- 
mally corresponds to E g E g configuration. In the present 
calculation there is no sign of any orbital ordering at 
all. We have calculated quadrupole moment tensor for 
d-shell of V ion and have found that while 3- fold symme- 
try is broken for "real" AFM structure, the quadrupole 
moment tensor is the same on all V ions. We should 
mention that we did not include the spin orbit coupling 
in these calculations and these could have appreciable ef- 
fects for V 3+ . For example the crystal structure allows 
for a Dzyaloshinsky-Moria coupling which could result 
in a small non colinear spin allignement. This could in 
the end increase the unit cell and could be visible in an 
annomalous scattering experiment. 

We should stress here that the total energy difference 
between different magnetic structures is really small. For 
the monoclinic crystal structure the "real" AF state has 
the lowest total energy followed by the simple AF state, 
which is higher by 80 K. In the case of the corundum crys- 
tal structure the simple AF state has the lowest energy, 
but the energy difference with the "real" AF state was 
only 5 K. Thus in the corundum phase these two mag- 
netic structures (real and simple) turned out to be almost 
degenerate. In short, the LDA+U results show us how 
important the small monoclinic distortions are for the 
magnetic structure and how low energy scale excitations, 
involving only spin reorientation, can effectively change 
large energy scale values such as band gaps. 

There are two problems arising from the present re- 
sults. First the spin should be considered to be 1 per 
V atom rather than 1/2 and secondly that the orbital 
occupation is consistent with a E g E g configuration of 
the ground state, which is in-plane symmetric and or- 
bitally nondegenerate, but still with the complex "Real" 
magnetic structure of the AFI phase. In this magnetic 
structure shown in Fig 1, every atom has three closest 
neighbors in the basal plane, one of which is ferromag- 
netically aligned (J a ) and the other two — antiferro- 
magnetically {Jp 1 and Jp 2 ), and also there is one neigh- 
bor along the hexagonal c axis which is ferromagneti- 
cally aligned (J 7 )(see Fig]l] for definitions of exchange 
interaction parameters). To check the above mentioned 
consistency we calculated the exchange interaction pa- 
rameters (EIP) using a well tested method described in 
p3| . We calculated EIP's for both crystal structures: 
monoclinic and corundum. It was found that according 
to EIP calculation only one stable magnetic structure ex- 
ists in the monoclinic phase. ("Stable" (or " consistent" ) 
means that if in the EIP calculations the certain pair of 
spin was parallel, the corresponding exchange parameter 
came out as ferromagnetic, and if antiparallel, as anti- 
ferromagnetic.) It is the "real" AF magnetic structure 
with the following values of EIP's: J a = 48 K (ferromag- 
netic), J Px = -214 K, J/3 2 = -90 K, J 1 = 47 K. For the 
corundum crystal structure we obtained two stable mag- 
netic configurations from the EIP calculation point of 



view: "real" AF (J a = 55, J Pl = J fh = -120, J 7 = 44) 
and magnetic structure with uniform AF exchange in 
the basal plane ( J a = Jp 1 = Jp 2 = —65) and small 
frustrated exchange along the hexagonal c axis (J 7 ~ 0). 
From these results we can say, that indeed the monoclinic 
distortion of the crystal structure stabilizes the real AF 
magnetic structure, and the E g E g configuration of the d 
electrons is consistent with this magnetic structure. The 
values of EIP's depend strongly on the magnetic struc- 
ture, which tells us that one cannot adequately model the 
magnetic interactions in V2O3 as only nearest neighbor 
Heisenberg exchange. This is most probably connected 
with the fact that V2O3 is close to being metallic which 
is also reflected in the strong dependence of the electronic 
structure on the magnetic one (and vice versa) which we 
saw in LDA+U calculation. The polarized neutron scat- 
tering experiments |p4| show the qualitative change of 
magnetic interactions in the transition from the antiferro- 
magnetic insulator with the monoclinic crystal structure 
to both the metallic phase and paramagnetic insulator 
with the corundum structure. Instead of a peak in re- 
ciprocal space corresponding to AFI magnetic structure 
("real" AF) it has a peak corresponding to a magnetic 
structure with all three V-V interactions in basal plane 
antiferromagnetic ("layered" AF). Another peculiarity of 
the neutron scattering results is a very large width of the 
peak. 

Our results show that for the distorted monoclinic crys- 
tal structure the "real" AF magnetic structure is the 
lowest one and the only one which gives a consistent 
set of exchange interaction parameters. Transition to 
the corundum crystal structure leads to the coexistence 
of the two well defined sets of exchange interaction pa- 
rameters with nearly degenerate values of total energy 
for those two magnetic configurations, the symmetric in 
basal plane "layered" AF structure being slightly lower in 
energy. That would result, in reciprocal space, in a peak 
centered at the corresponding Q- vector but this peak will 
be strongly broadened due to the transitions to the ex- 
cited state and corresponding fluctuations in the mag- 
netic structure. 

Our calculations also give a value of the magnetic mo- 
ment per V ~ \.l \XBi nearly the same in all the structures 
studied. This value is somewhat larger than the value of 
1.2 [lb obtained from neutron scattering for the antiferro- 
magnetic phase, but is consistent with the value (1.7 \ib) 
obtained from the high-temperature susceptibility. This 
relatively large value of \i is evidently a consequence of 
a strong Coulomb interaction on V 3+ which tends to de- 
stroy the formation of the molecular orbital singlet state 
on j4i g -orbitals of V-V pair, which was assumed in most 
previous studies. The fact that this is independent of the 
magnetic or crystal structure suggests that this should be 
treated as a high energy scale parameter in any model. 

Summarizing, we have carried out LDA+U calcula- 
tions of the electronic structure and exchange constants 
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of V2O3 in both the monoclinic and corundum structures 
and obtained a consistent description of the main proper- 
ties of the antiferromagnetic insulating phase and of the 
paramagnetic insulating one. In contrast to the previous 
assumptions, in both these phases the electronic config- 
uration is predominantly E g E g one, i.e. two d-electrons 
of V 3+ occupy the doubly degenerate _E 9 -orbitals. In ad- 
dition the spins of the two electrons are parallel leading 
to a high spin S=l local moment. As a result there is 
no orbital degeneracy left and correspondingly no orbital 
ordering of the kind which was invoked previously to ex- 
plain magnetic properties of V2O3 [^10[. Despite that 
we are able to obtain the correct magnetic structure of 
V2O3: the signs of the exchange constants in the mon- 
oclinic phase are consistent with the observed antiferro- 
magnetic structure. The calculated values of the energy 
gap ~ 0.6 eV and the magnitude of the magnetic moment 
per V ~ 1.7jU_b also agree with the experiment. Strong 
change of magnetic correlations through TV is ascribed 
to the near degeneracy of different magnetic structures in 
the corundum magnetic structure of V2O3, the state with 
antiferromagnetic correlations in all three directions in a 
basal plane having slightly lower energy. Thus orbital or- 
dering is not required to explain the physical properties of 
V2O3 in the antiferromagnetic and in the paramagnetic 
insulating phases. The LDA+U results also demonstrate 
that a spin 1/2 Hubbard model is not the correct start- 
ing point but that one should use a spin 1 model with 
very strong hunds rule exchange. Within such a model 
with two electrons S=l in E g orbitals on each site the 
hopping would have to involve both the minority spin 
Eg orbitals as well as the A% g orbitals which are close 
in energy. Having found this it is also not so surprising 
any more that the electronic structure and especially the 
d band widths and splittings are strongly dependent on 
the nearest neighbor spin spin correlation functions and 
that there may be a redistribution of electrons between 
the Eg and Ai g states on going into the PI or Metallic 
phase as found in recent electron spectroscopic studies 
& _ 
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FIG. 3. Partial densities of states for d states of V ob- 
tained in LDA+U calculation, a) "Real" antiferromagnetic 
structure; b) "simple" antiferromagnetic structure; c:) ferro- 
magnetic structure 
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